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(57) Abstract 

The invention relates to a receiving method and a receiver in a system in which received signals comprise a number of summed 
signal components of which at least some are unknown, and the receiver comprising a means (202) for sampling the received signal and 
forming an information vector of the samples, and a means (208) for estimating the received signals with one or more multidimensional 
filters. To enable reliable and advantageous signal detection, the receiver comprises a means (208) for forming a model of the signal by 
means of a cost function which reveals non-linearities in the signal, and for optimizing a filter function, and a means (208) for detecting 
the data by means of the optimized filter functions so that prior to the detection the dimension of the information vector is reduced. 
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Receiving method and a receiver 

FIELD OF THE INVENTION 

The invention relates to a receiving method in which a received 
signal comprises several summed signal components, at least some of which 
5 being unknown, and the received signal being sampled, said samples forming 
an information vector, and in which the received signals are estimated with 
one or more multidimensional filters. 

PRIOR ART 

In planning and implementing communication systems, one of the 

10 fundamental problems is concurrent transmission and reception of signals by 
several simultaneous users so that the signals would interfere with one 
another as little as possible. For the above reason, and because of other 
interference caused by the transmission channel, as well as the transmission 
capacity being employed, a number of different transmission protocols and 

15 multiple access methods have been developed. Particularly in mobile 
communication systems, the most common of these include the FDMA and 
TDMA methods, and lately also the CDMA method. 

In all the above methods, interference is present in somewhat 
different forms, and different kinds of methods have been employed in these 

20 methods to overcome interference-related problems. In TDMA, a significant 
source of interference is intersymbol interference (ISI), and in CDMA, multiple 
access interference (MAI). Generally speaking, in detection of a desired 
signal, the aim is to find the desired symbols in a signal to which surplus 
information has been summed. The different multiple access methods pursue 

25 this aim in different ways. For CDMA, simultaneous detection of a plurality of 
users, and for TDMA different types of equalizers have been suggested. 

The equalizers, which seek to alleviate signal distribution in time 
domain, come in two basic types: linear and non-linear. The linear equalizers, 
such as Wiener filters, are simple to implement but do not perform well enough 

30 in all situations due to a modulated waveform being characteristically non- 
linear. Non-linear equalizers are represented by e.g. a Viterbi equalizer and 
decision-feedback equalizers. 

Correspondingly in the CDMA, the MUD methods are either linear 
or non-linear. Of the former, a decorrelating receiver, and of the latter, a Viterbi 

35 and decision feedback coupling method may be mentioned. A closer 
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description of decorrelating receivers may be found in e.g. A. Hottinen: Self- 
organizing multiuser detection, Proc. of IEEE ISSSTA '94, Oulu, Finland. 

Also known in prior art is the use of so-called blind adaptive 
algorithms in the equalizers and detection. A blind algorithm refers to a 
5 adaptation rule that does not take into account the polarities of the signals, i.e. 
a training sequence is not required to bring the receiver in a correct state. 
Blind methods are described in association with CDMA, for example in the 
reference S. Verdu: Adaptive multiuser detection, Proc. of IEEE ISSSTA *94, 
Oulu, Finland. 

10 All the above methods, including the non-linear ones, share the 

drawback that as the number of interfering signals increases, which means 
that the dimension of the computation space increases, a steep increase 
occurs in the computation amount required for carrying out optimal detection. 
As the optimal detection has thus turned out to be difficult, suboptimal 

15 methods have been employed, which has resulted in loss of received 
information, manifest at the detector as deteriorated signal quality. 

CHARACTERISTICS OF THE INVENTION 

It is consequently an object of the present invention to implement a 
method by means of which a received signal containing interference may be 
20 detected suboptimally even in case the signal model is unknown. The 
invention further aims at implementing an adaptive method which requires little 
information on the interfering signals and noise and which is able to detect the 
desired signal reliably with a smaller amount of information even in a disturbed 
propagation environment. 
25 These objects are achieved by a method of the type set forth in the 

introduction, characterized in that signal filtering takes place adaptively by 
forming a model for the signal by means of a cost function which reveals non- 
linearities in the signal and which is used in optimizing a filter function, and 
that the data are detected by means of the optimized filter functions so that 

30 prior to the detection the dimension of the information vector is reduced. 

The invention further relates to a receiver in a system in which 
received signals comprise a number of summed signal components of which 
at least some are unknown, and the receiver comprising a means for sampling 
the received signal and forming an information vector of the samples, and a 

35 means for estimating the received signals with one or more multidimensional 
filters. The inventive receiver is characterized by comprising a means for 
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forming a model of the signal by means of a cost function which reveals non- 
linearities in the signal and for optimizing the filter function, and a means for 
detecting the data by means of the optimized filter functions so that prior to the 
detection the dimension of the information vector is reduced. 
5 The method and receiver according to the invention provide a 

multitude of advantages. The receiving method of the invention employs so- 
called projection pursuit, and it is capable of efficiently removing the effects of 
intersymbol interference and multiple access interference, and in general 
nonlinear interference in the signal to be detected. It is applicable in 

10 connection with FDMA, TDMA as well as CDMA. It may further be used in a 
multi-stage receiver as one of the latter stages. For example, in connection 
with the CDMA, a conventional receiver stage implemented by filters adapted 
to spreading codes may serve as the first receiver stage the output signal of 
which is processed by the receiver stage according to the invention. A large 

15 amount of computation associated with the prior art methods may in the 
solution of the invention be avoided by restricting estimation to some of 
subspaces chosen, without, however, losing any significant information. 
Further still, the method makes it possible not to use a correlator in the multi- 
stage detector for the interfering users. 

20 As mentioned, the solution according to the method of the invention 

utilizes a so-called projection pursuit method which may be used for finding 
the ideal filters required for signal processing. The projection pursuit scheme 
applies a general projection index, i.e. a cost function, used to find a subspace 
indicating the most significant features of a signal. According to the method, 

25 the receiver further finds a number of subspaces that together collect 
information on a time or code correlated signal, this information being utilized 
in the detection. The projection pursuit method is known per se in statistics, 
and it is described e.g. in P. Huber: Projection Pursuit, Annals of Statistics, 
Vol. 13, No. 2, pp. 435-475, 1985, which is herein incorporated as reference. It 

30 should be noted that the solution according to the invention applies projection 
pursuit in a new way in a telecommunications context in which the method is 
used to detect the nonlinear features of a signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following, the invention will be described in closer detail with 
35 reference to the examples in the accompanying drawings, in which 
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Figure 1 illustrates a cellular radio system in which the method 
according to the invention may be applied, 

Figure 2 is an exemplary illustration for the structure of a receiver 
implementing one of the embodiments of the invention, 
5 Figure 3 illustrates projecting a set of BPSK signals on an x-axis, 

and 

Figure 4 illustrates an example of implementing projection pursuit. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 is an exemplary illustration of a typical cellular radio 
10 system in which the receiving method according to the invention may be 
applied. In this exemplary case, a CDMA system is employed but the invention 
is applicable to the other multiple access methods, as is obvious for a person 
skilled in the art. 

The system comprises a base station 100 and a number of 
15 subscriber terminals 102 - 106 which have a connection 108 - 112 to the base 
station. In case of a CDMA system, all the subscriber terminals 102 - 106 
communicate with the base station 100 on the same frequency band. The 
transmissions from the terminal equipments are distinguished by a pseudo- 
random spreading code. Although the aim is to choose the spreading codes 
20 used on the different connections so that in an ideal case they do not correlate 
with each other, i.e. do not interfere with each other, in practice the other 
connections are, during reception of the desired signal, perceived as multiple 
access interference. The interference may originate either from the same or 
surrounding cells. The method of the invention is applicable to the receiver of 
25 both the base station and the terminal equipment. 

In the following, an exemplary description is at first offered of an 
asynchronous CDMA system in which the method of the invention may be 
applied. As is obvious for a person skilled in the art, the invention may also be 
applied to a synchronous system. 
30 u is assumed that there are K users in an asynchronous BPSK 

CDMA system, each having been assigned the spreading waveform 

su(t) = tsmr.O-jTj 
j- 1 

in which the /th chip, i.e. bit, of the spreading waveform for the Mh 

user is obtained from 

35 St e {-U},k = l K. 
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In the above, T c is the chip length and N signifies the number of 
chips in the spreading waveform. No restrictions are set for a user's spreading 
code. The users transmit the information by modulating the spreading 
5 waveforms with the data symbols bff e A in which A is a symbol alphabet, in 
the case of BPSK {-1,1}. In CDMA, then, all the users transmit simultaneously 
on the same frequency. It is therefore the task of the receiver to demodulate 
the summed signal 

k p Ek 

r(0 = Z H h "0)*H\— s k (t-i1) + n(t) 
k-i i-p | T 

10 

where n(t) represents white Gaussian noise at which the power spectral 
density HJ2, 2P+1 denotes the length of the data packet, h„(t) denotes the 
impulse response of the Wh user's physical channel, E k denotes the 
transmission power of the tth user. The symbol * denotes convolution. The 
15 signal propagates through a multipath channel, and the impulse response h^t) 
of the channel may be represented by the formula 

i 

1*0 

20 where the fth channel coefficient, i.e. the channel tap, of the *th user is h kt x e 
c and t is the corresponding delay. When the coefficient is time-varying with 
Rayleigh distributed amplitudes, the channel is called a Rayleigh channel. 
When h K , is constant, the channel is referred to as a Gaussian channel 
(assuming that the noise n(t) is Gaussian noise). 

25 The receiver receives the signal r(t), and correlates the signal with a 

receiver filter w k (l), i.e. 

(l+l)T*t u 

z% = J r(t)wUt-T U -iWdt 

IT*t u 

30 where T denotes the length of the symbol interval. In CDMA, the filter is 
typically matched to the transmitted waveform, w k = s k . 
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In order to build a matrix model of the received signal, the outputs 
of the Ath user's matched filters are collected to the vector 

if - (& &) 

5 the outputs of all K users' filter outputs to 

z«> = (zf m 

and finally a sequence of matched filter outputs to 

10 

z = (z f ....,z r ? 

The amplitude related to the Ah symbol of the Mh user's Ah path is denoted 
with a$ . The amplitudes 

15 ccv(0 s \Y 

are assumed to be constant during a symbol interval, and 

20 

The amplitudes of the Mh user's channel coefficients at the Ah time interval are 
collected to 

of = r<«* • 

25 The matrix of amplitudes for signals received by all K users is determined as 

A°> = dag (off a ( D 

30 The correlation between the Mh user's Ah multipath component and the fclh 
user's /th multipath component is obtained from 

CO 

[RSl-Jtr = f w k (t-Tu)Sf(t-ttjr + nT)d 
By collecting these to matrices 
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5 Rs = 



R% - R ( 2k 
r«» ... J?™' 

• • • 



the received signal may be written in the form 
10 z = RAb + n, 

where R ■ Rn, a = diag (A ( - 2N> , . . . ,a ,2n) ) and b = (b 

<- E \...,b ,p, ) T , where b (i > = (bfl T m T ) . Here, l = (l,...,l) T is 

an L vector. The noise is assumed Gaussian with covariance R. 
15 The parameters of the received signal, i.e. the delays and 

amplitudes, are estimated in the receiver, and assuming that the signal has 
propagated in a Gaussian channel, the prior art parametric detectors, such as 
ML (Maximum Likelihood), solve 

20 b op t = arg min \z - Rab| 2 c , 

where |y| 2 c = y M Cy ja c is the weighting matrix. Here, c = (R+8 2 A)" 1 . 

According to the formulas above, a prior art receiver detects the 
desired signal from among interfering signals, providing the computation load 
25 does not increase too much or the signal model is known. The method 
according to the invention is particularly well applicable in case the model is 
known in part only (a submatrix of R) or if the noise is not Gaussian 
distributed. 

In the following, the adaptive method according to the invention will 
30 be described. It should be noted that almost all receivers are adaptive to some 
extent. Different adaptive features are determined by which information the 
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receiver assumes to know of the signal a priori, what it has to learn about the 
received signal, and what it will do with the information learned. The solution 
according to the invention is based on the statistical decision theory having its 
foundation on a criterion that is dependent on the probability density function 
of the signal. If the receiver knows the density function, it is able to make an 
optimal decision on the basis of e.g. ML, MAP (Maximum A Posteriori) or a 
general Bayesian criterion. A key feature of the invention is that the probability 
density functions are estimated directly from the received signal. 

In the projection pursuit method, the estimated signal is processed 
with one or more multidimensional filters, and by means of the cost function 
the receiver adaptively finds filter functions for the received signal, the cost 
function distinguishing nonlinear features, such as non-Gaussian features or 
clusters, in the signal. The receiver typically consists of several stages, each 
stage being able to carry out dimension reduction and signal detection. The 
15 detector may be implemented by various ways, for example by nonparametric 
density estimation. The nonparametric method may efficiently be implemented 
by utilizing Fourier transformation (FFT). 

Figure 2 is a block diagram illustration of a simplified structure for 
the front end of a receiver implementing one of the embodiments of the 
20 invention. By means of a radio frequency means (not shown), transducing 
means 202 are operationally connected to an antenna 200. The received 
signal 216 is applied to a converter 202 which converts the signal to digital 
form by sampling it at a desired sampling frequency. The sampled signal 210 
is applied to a memory element 204. The samples are read out from the 
25 memory element 204, and are fed to a signal processing processor 206 which 
performs a linear transformation to the signal, e.g. permutation. Interleaving, 
preliminary dimension reduction, or decoding. The permutation may be the 
inverse function of the permutation performed e.g. at the transmitting stage. 
Samples may be taken out from the processing means 206 at a desired rate. 
30 Dimension reduction may be performed so that the output rate of the 
processing means is lower than the input rate. However, the memory element 
204 and the processing means 206 are not mandatory, and the solution 
according to the invention may also be implemented without them. 

From the processing means 206, the signal 214 is applied to a 
35 projection pursuit means 208 which processes, detects, and estimates the 
constituent symbols either parametrically or nonparametrically, or their 
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combination. Nonparametric estimation refers to a method in which no initial 
assumption is made of the signal model or the interference source, or of the 
density estimate. The latter way is inefficient in high dimensions, and therefore 
subspace estimation is used for dimension reduction. As a person skilled in 

5 the art will know, the inventive receiver also comprises other components, 
such as filters. For reasons of clarity, these components are not dealt with 
here as they are inessential from the point of view of the invention. 

Assuming that { P m } , m= l , . . . , m is a set of projection operators 
to subspaces {L m e c*"} , m=l,...,M, d-l,...,D and d denotes 

10 dimension. Here, the subspaces may be of different dimensionality, and may 
be implemented by a set of 2d ro filters that may be adaptive of fixed. The 
method may be described with the aid of the following steps: 

■ Find a set of low dimensional subspaces from the received signal given, by 
using suitable criteria or cost function. 

15 ■ Estimate statistics of signal r, in the low dimensional space (density or the 
parameters of the density). 

■ Remove the information in the low dimensional space (by using previously 
estimated functions or/and parameters for transformation of signal) to obtain 
processed signal r M . 

20 ■ Combine the statistics from the previous dimensions and receiver stages. 

■ Go back to step 2. 

What is essential, then, is to find the low dimensional subspace L. 
The cost function used in the projection pursuit method of the invention 
reveals nonlinear features in the information. Typical non-linear features (non- 
25 Gaussian signals) are formed e.g. in the modulated data. To take an example, 
when the set of BPSK signals is projected to the x-axis, a clustered data set is 
formed because it is a combination of two Gaussian distributed distributions 
with different means. This is illustrated by Figure 3. The same principle applies 
in higher dimensions as well. In the communication context, higher dimensions 
30 are formed when the data is spread e.g. in time according to channel impulse 
response, or in symbol domain (like in CDMA). 

Projection pursuit provides means for collecting information from 
the multidimensional signal by finding characteristic features of the signal and 
by combining them effectively. The features are revealed by filters optimized 
35 adaptively by using cost functions different from the ones typically used in 
communications contexts. 
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The cost function should be defined so as to reveal the interesting 
features of the signal. There are two basic ways to choose the cost function: 
the statistical features of the density (e.g. moments, cumulants) or the density 
estimate as such. The density, in turn, may be estimated from the statistics or 

5 functions of the statistics by parametric density estimation, or directly from the 
samples by non-parametric density estimation. 

Many distance measures between densities may be employed as 
cost functions. One of them is the so-called F-divergence. F-divergence may 
be used as the cost function in an adaptive projection pursuit receiver. 

! o For example, if f(u) = ulog(u), then 



F+* = IP.wfjdA 

where p, and p 2 are densities, whereby the F-divergence is the same as the 
15 Kullback-Leibner distance. 

To take another example, if f(u) = 1 1 -u 1/p I p , then 

20 whereby the F-divergence coincides with Hellinger distance. 

Other measures may also be used as cost function, depending on 
the features to be revealed from the signal. Such other measures are 
described in the aforementioned publication by P. Huber. 

If particularly the CDMA system and multiuser detection are 
25 examined, the non-Gaussian features of the signal can be seen e.g. in the 
residual signal from which the effect of known interference has been 
cancelled. It is possible to estimate the effects of unknown interference by 
means of a suitable cost function on the basis of the nonlinearity caused by 
the interference. Hence, interference cancellation may be carried out without 
30 having to detect the interfering bits. 

In the following, the projection pursuit method according to the 
invention is examined with the aid of another example. The desired subspace 
is assumed to be and the corresponding projection matrix P^. It is further 
assumed that M samples have been stored in the matrix z = 
35 ( z 111 , . . . , z m ) ■ m denotes the so-called kernel, i.e. the probability density 
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function. In order to determine the subspace. it is first necessary to determine 
the cost function Q(.) that finds the desired features. In case of nonlinearities, 
the cost function Q should be invariant to afRne transformations, Qi (AZ+m) = 
Qi (Z) .An affine transformation refers to the combined result of one or more 
5 parallel projections. This cost function class is able to capture higher order 
features of the data. Cost functions that satisfy q 2 (AZ+m) = aq 2 (Z) +m are 
suitable for extracting first order information (mean, location), and the ones 
that satisfy Q 3 (AZ+m) = AQ 3 (Z)are sensitive to second order information 
(covariance). 

10 The following exemplary algorithm illustrates the use of projection 

pursuit. Denote z(0)aZ and Q(.)= Qi - f is an unknown d-dimensional 
probability density function, 

(1) Set Z(i+1)- Whiten (Z(i)). 

(2) Solve P^O) = argmax p Q(PZ(i+1)). 
15 (3)SetY(M) = P w- £(i+1). 

(4) Update the model f l+ i = f ih i+ i (Y (i+l ) ) . 

(5) Generate a sample from the new model z(i+l) - fm. 

(6) Go to step 1 . 

20 Step 1 is optional and does not affect the optimization for the 

selected cost function class. The whitening provides the advantage that all 
marginal densities will have variance 1, and only the higher order information 
remains in the data. This is useful in implementing step 2 for computational 
reasons. Steps 1 and 5 may be implemented in many ways. For example, a 

25 Monte Carlo sample may be taken from the new model as in analytical 
projection pursuit, or the new sample may be taken by using a non-linear 
transformation. The interference cancellation step in the multistage detector 
provides one non-linear method for Gaussinizing the data. The residual signal 
z-Rab(m) of the interference cancellation is Gaussian with zero mean if the 

30 interference cancellation succeeded. 

The projection pursuit scheme may be used for finding, from 
multidimensional data, such a projection level at which the signal is non- 
Gaussian, and in the subspace in question the data may be Gaussinized by a 
decision-feedback coupling, as in interference cancellation, by the Monte 

35 Carlo method or non-linear transformation. This is explained e.g. in J. 
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Friedman: Exploratory projection pursuit, J. American Statistical Association, 
Vol. 82, No. 397, pp. 249-266, March 1987. 

In the following, the implementation of the algorithm described 
above is examined in the receiver by means of the block diagram in Figure 4. 
5 The figure shows how a signal is applied, in the first stage 400 of the 
projection pursuit, first to an adaptive filter 402, which may be implemented 
e.g. on the rake principle. From the filter, the signal is applied to a processor 
404, in which signal whitening is performed. The whitened signal is applied to 
a processor 406, which finds non-linear features, i.e. computes P^. Next, the 
10 density estimate is updated in a processor 408 and a sample is generated 
from the new model in a processor 410. The sample z(1) thus obtained is 
applied to the second stage 402, is which the procedure is repeated, and a 
new sample z(2) is obtained, etc. The aforementioned processors may 
naturally be replaced in practice with one signal processor into which the 
15 algorithm required has been stored. As far as the invention is concerned, the 
whitening mentioned above in block 404 is not essential. The estimate for the 
density function is corrected on the basis of data z(0), z(1),.., z(n). 

As known, detection generally relates to recognizing user bits 
among external interference in the received signal. The detection is based on 

20 the signal z and the estimated density from successive sequences z, z«. 

After the densities have been estimated, the detection employs optimal 
Bayesian rule. 

Although the invention has in the above been descnbed with 
reference to the examples in the accompanying drawings, it is obvious that the 
25 invention is not restricted thereto but may be varied in many ways within the 
inventive idea of the attached claims. 
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CLAIMS 

1 . A receiving method in which a received signal comprises several 
summed signal components, at least some of which being unknown, and the 
received signal being sampled, said samples forming an information vector, 

5 and in which the received signals are estimated with one or more 

multidimensional filters (208), characterized in that 

signal filtering takes place adaptively by forming a model for the 

signal by means of a cost function which reveals non-linearities in the signal 

and which is used in optimizing a filter function, and that 
10 the data are detected by means of the optimized filter functions so 

that prior to the detection the dimension of the information vector is reduced. 

2. A method as claimed in claim 1, characterized in that the 
data transfer employs code division multiple access, and that the initial values 
of a multidimensional filter are computed by means of filters matched to 

15 spreading codes used on each connection. 

3. A method as claimed in claim 1 or 2, characterized by 
performing, prior to filtering, parametric detection for the known part of the 
signal, and by subtracting this information from the received signal. 

4. A method as claimed in claim 3, characterized in that the 
20 parametric detection comprises decorrelating detection. 

5. A method as claimed in claim 1, characterized in that the 
cost function reveals non-Gaussian features in the signal. 

6. A method as claimed in claim 1, characterized in that the 
signal model is employed for estimating the density function of at least one 

25 signal interfering with the desired signal, and that the estimation of the density 
function takes place at subspaces determined by the filter functions, and that 
the detector combines the density functions of the subspaces. 

7. A method as claimed in claim 1 or 2, characterized by 
carrying out a linear transformation for the signal prior to the multidimensional 

30 filtering. 

8. A method as claimed in claim 7, characterized in that the 
linear transformation comprises the steps of storing signal samples in a 
memory element (204), and reading out the samples from the memory 
element in a different order than in which they were stored. 

35 9. A method as claimed in claim 1, characterized by 

employing the cost function for finding, for the multidimensional data, a 
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subspace in which the data are non-Gaussian, and by Gaussinizing the data 

distribution in said subspace. 

10. A receiver in a system in which received signals comprise a 
number of summed signal components of which at least some are unknown. 
5 and the receiver comprising a means (202) for sampling the received signal 
and forming an information vector of the samples, and a means (208) for 
estimating the received signals with one or more multidimensional filters, 
characterized by comprising 

a means (208) for forming a model of the received signal by means 
10 of a cost function which reveals non-linearities in the signal, and for optimizing 

the filter function, and 

a means (208) for detecting the data by means of the optimized 
filter functions so that prior to the detection the dimension of the information 
vector is reduced. 

! 5 1 1 A receiver as claimed in claim 1 0, c h a r a c t e r i z e d in that 

the receiver comprises a means (402) for carrying out adaptive filtering for the 
signal a means (404) for whitening the filtered signal, a means (406) for 
revealing nonlinearities in the signal, a means (408) for computing and 
updating the density estimate, and a means (410) for generating a new 

20 sample by means of the updated density estimate. 

12. A receiver as claimed in claim 10, c h a r a c t e r i z e d in that 
the receiver comprises means (204, 206) for carrying out a linear 
transformation for the sampled signal (210) prior to the multidimensional 

filtering (208). m 
25 1 3. A receiver as claimed in claim 10, characterized in that 

the means (208) comprise one or more estimating stages (400, 402). 
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